Single or multi-flavor Kondo effect in graphene 
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Abstract 

Based on the tight-binding formalism, we investigate the Anderson and the Kondo model for an 
adaom magnetic impurity above graphene. Different impurity positions are analyzed. Employing 
a partial wave representation we study the nature of the coupling between the impurity and the 
conducting electrons. The components from the two Dirac points are mixed while interacting with 
the impurity. Two configurations are considered explicitly: the adatom is above one atom (ADA), 
the other case is the adatom above the center the honeycomb (ADC). For ADA the impurity is 
coupled with one flavor for both A and B sublattice and both Dirac points. For ADC the impurity 
couples with multi-flavor states for a spinor state of the impurity. We show, explicitly for a 3d 
magnetic atom, d z 2, (d xz ,d yz ), and (d x 2_ y 2,d xy ) couple respectively with the Ti, T^(Ei), and Tq(E2) 
representations (reps) of Cq v group in ADC case. The basses for these reps of graphene are also 
derived explicitly. For ADA we calculate the Kondo temperature. 

PACS numbers: 75.20.Hr,72.15.Qm,71.55.-i,81.05.ue 
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INTRODUCTION 

In Graphene, a monolayer of carbon atoms recently fabricated successfully []] , the valence 
and the conduction bands touch at two inequivalent Dirac points K_ and K + at the corners 
of the first Brillouin zone (FBZ). Near K_ and K + the low energy dispersion is linear, 
indicating a massless Dirac fermions behavior. This particular band structure is at the 
heart of a number of unusual electronic properties [2j. Graphene is also an interesting 
candidate for transport-applications, in particular for spintronics: It exhibits a remarkably 
high mobility and the carrier density is controllable by a gate voltage; the mean free path 
can be as large as 1 fim. Graphene is however not immune to disorder that influences 
its electronic properties [2j. Extrinsic disorder is realized in a variety of ways: adatoms, 
vacancies, charges on top of graphene or in the substrates, and extended defects such as 



cracks and edges. When localized magnetic impurities are added [3| the Kondo effect, i.e. 
the dynamic screening of the localized moment, has to be addressed at temperatures T 
below the Kondo temperature T K . In this context, previous studies addressed the influence 
of magnetic impurities using the Hatree-Fock approximation 0, 0] which is valid at T > 
(see also fl). In Ref. [J the anisotropic single channel Kondo model is investigated briefly 
and in Ref. jsj the inftnite-U Anderson model for an impurity embedded in a graphene 
sheet has been employed and concluded a Fermi liquid behaviour. Recently, it has been 
claimed that a two-channel Kondo in graphene is present due to the valley degeneracy of 
the Dirac electrons [3], which leads to an over-screening and thus to a non- Fermi-liquid- 



like ground state. Very recently Refs. lCj, LL3J reported on scanning tunneling spectroscopy 



studies to investigate the Kondo effect in graphene. All of these studies consider a particular 
configuration of the impurity. In this work, we show that the position of the impurity on or 
in graphene plays a subtle role and affects essentially the underlying physics. In addition to 
the role of the vanishing density of state (DOS) at the Dirac points and the linear spectrum 
in their vicinity, a further important issue is to clarify whether a single or a multi-channel 
Kondo problem is realized. This is insofar important as the ground states for these two cases 
are essentially different: As established [121 ]. when the channel number iV cha is equal to 2S ira 
(£i m is the spin of the impurity), the impurity spin is then compensated by the host electrons 
completely, resulting in a singlet Fermi liquid ground state. When > 2S' im we enter the 
over-compensated regime where an opposite spin to the original spin appears and acts as 
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(a) (b) 

FIG. 1: (color online), (a) shows graphene primitive vectors and the symmetry operations defining 
the group of the crystal, a v i,i = 1,2,3 are reflection planes showing up in Cq v little group at T 
point, (b) shows FBZ with the symmetry operations (see text for explanations). 

a remainder spin; an effective antiferromagnetic coupling between this remainder effective 



impurity spin and the conducting spins results in a non- Fermi liquid ground state [131-115]. 
Thus, it is crucial to clarify when graphene with an impurity with a spin one-halve is a 
single or a two channels Kondo system. Recently, scanning tunneling conductance spectra 
are investigated with respect to the position of the scanning tip in Refs. [if], Q (see also 
Ref. |ia|). 

As a first essential step we identify the physical origin of the multi channels for the Kondo 
effect by analyzing the tight-binding Anderson model 19|, |20j. A partial wave method is 
then employed to reduce the redundant flavors and to identify the coupling flavors. The role 
of the Dirac points is exposed. An adatom above one carbon atom (ADA) or an adatom 
above the center of the honeycomb (ADC) configurations are considered. 

THEORETICAL FORMULATION 

We start from the Anderson Hamiltonian 

H = H g + Hf + H hyh . (1) 
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H g (Hf) describes graphene, (impurity) and i^hyb stands for the impurity-graphene hy- 

— 2 )' Wnere 



■1 Vz\ 



> 2' 2 a 2 a V2' 2 



Voa cc , a cc is the distance of the nearest carbon atoms. The B sublattice is related to 



bridization. The graphene primitive vectors are ai 
a 

A sublattice by Si = (a x — a 2 )/3, S 2 = ax/3 + 2a 2 /3, and S s = —Si — S 2 = — 2a x /3 — a 2 /3 
(cf. Fig. 1(a)). The two inequivalent Dirac points in FBZ are K.± = ±—(2/3, 0). The second 
quantization tight-binding-Hamiltonian 0, [21 23 1 reads H g = Y^imafalabja + t*b^ a a ia ) 
where the sum runs over the nearest-neighbor pairs (ij), o is a spin index, and ai a (bi a ) 
are annihilation operators for states on the A (B) sublattice. The Hamiltonian in momen- 
turn space is H g = £ fc fc(k)ak&fc + r(k)6^o E(T ], where £(k) = -t$(k), $(k) = £ 5j e** 
and k are wave vectors in FBZ. The hybridization Hamiltonian in tight-binding formalism is 
#hyb = E^[(Vfalf a + V^blf a )+h.c), where = -fe J2^ m e-*<*^V«t(Jl t + 
T a — Rim), a=A or B, T a stands for the relative location for a atom in the unit cell, 
V a f(Ri + T a — Ri m ) = f <ir[0 a (r — Rj — r Q )]*/i(r)0i(r — R im ), and (fi a is the atomic function 
for the atoms in a sublattice (only 7r orbitals are included), <pi is the localized impurity 
wave function, h is the single particle Hamiltonian, and N is the number of unit cells. 

V a f(Ri + T a — R; m ) is a Slater-type bond in LCAO formalism 24j which describes 
the strength of hybridization between two atomic orbitals located at the impurity and its 
neighbors. Here only p n electrons are relevant. They form the ir and ir* which touch at the 
Dirac points K±. An s-wave impurity substituting one carbon atom in the graphene plane 
is decoupled from tc and 71* bands; for d- and f-wave impurity in this substitution case is 
possible. Our focus is however on the case where the impurity is above the graphene plane. 
The relative positions of the impurity determine the hybridization strength, i.e. V a ^ and 
the phase is included in the exponential in V? . 

To obtain an effective low-energy Hamiltonian we expand the wave vector k = K-t + k 
around K ± . Thus, $(k) ps (—v F /t)(±k x - ik y ), or $(k)| s « (-v F /t)k Xe' iXe , where A = 
sgn(s) = 1(— 1) for s = K + (K_), v F is graphene Fermi's velocity, 9 is the azimuthal angle 
of k. The Hamiltonian of graphene is then 
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H g is diagonalized by introducing the ( fields [25] as 

« skCT = ^[C;.(k) + C CT (k)], 6 sk . = ^e^[C; CT (k) - C^k)], 
H g = E [^l fc lC;t(k)C; CT (k) - ^iCUk)C CT (k) 

scrk 



(3) 



Partial wave method and hybridization Hamiltonian 



Going over from 



discrete to 



£ SCT f dk ^|fc|C(k)C; CT (k) -^|A;|Ci(k)C CT (k) 



continuum 



26] 



we write 



and ifi 



hyb 
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Qo is the unit cell area. In 2D 



orbital momentum eigenfunctions e imd we write C± CT (k) 
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where {(™ s (k), (™ a s , \k')} = 5 ss >5 vv i5 aa >5 mm >5(k - k'). H g is then expressed as 

POO _. 

H g = j2 v F \k\dk c:\k)c™(k)-e: j (k)c:(k) . 



(4) 

For discussing the hybridization Hamiltonian for various geometric configurations we set the 
spatial zero point "O" as the projection of impurity position onto the unit cell. 
Adatom impurity above one A atom (ADA): In this case ta = 0, in the Dirac cones we find 
(we assume vq is real) V^f = V Jj^^ = The next nearest neighbors are three B atoms 



at <M = 1,2,3. Tbm,V*' = Zy£L 
we find 

-f^hyb 



A 



'■^= , where v% = V s $ . Substituting in H) 



hyb 



E 



'nfto\k\dk 
2^ 



vo(C s Ak)+Cf(k)) + 



fa + h.C. 



This equation tells us only m = flavor couples to the impurity. 

Adatom impurity above the center of the honeycomb: Here, the impurity hybridizes with the 
same strength t> , with 6 nearest neighbors, 3 A atoms and 3 B atoms. With respect to O, 
A2, B2, A3, B3 are in other unit cells. Three A(B) atoms yield $*($) when the phases are 
coherently superimposed, i.e. V^a^fa+V^ b\. a f a = 



VQ 



e -ik-5 3 (l _|_ e -ik.(a 1+ a 2 ) + g-ik-aij a ^J^ + g <k.*i ^ + g-ik-ai + e ik-a 2 ) 



$*(k)aL/, + $(k)&L/-] • For ^hyb we deduce 
wov^, f°° ^\k\dk 
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(5) 



When s = K + , A = 1, the impurity couples with the partial wave m = 1 from A sublattice, 
and m = —2 from B sublattice. When s = K_,\ = — 1, it interacts with m — — 1 from A 
sublattice, and m = 2 from B sublattice. 



FLAVOR RIGHT MOVERS AND SYMMETRY ANALYSIS 

We unfold the range of momenta k from (0, oo) to (— oo, +oo) by defining fla- 
vor right movers. For ADA case applies c\ a (k) = C+o-G^Dj f° r k > and 
clM = C°4(\k\), for A; < 0. We obtain thus H g = J2 sa I^L £ kdkc 1 ] 7 (k)cl (7 (k) : H hyh = 



couples wit 



where e k = hvpk. So, the impurity 



Too 4^ (-0 + ) ct{k)U + h.C. 

;h only one flavor of conducting waves. For ADC case the impurity couples 
with m = A or m = — 2A flavors from A or B sublattices respectively. Therefore, we 
introduce two flavor right movers as cf CT = \(+*(\k\), c s 2(J =C+l X ' s {\k\) for k > 0; c\ a =- 
ACi?(l*l), 4,=C^(l*l) for < 0. In terms of these, # g =En S . IZ e k dkc s J a (k)c s n(T (k), and 



LI ^^c* (*)/, + h.c. 
have {^(A;),^,^')} = 8 nn i5 ss i5 aa i5(k - k'). 



^hyb- ( _ t) 2^ 



-oo 

, where n = 1,2. In these two cases, we 



Symmetry analysis 

The valence bands formed by a bonds of sp 2 hybridization were discussed by Lomer long 
ago . The popular tight-binding Hamiltonian for ir and n* bands formed by p z orbital 



located on each carbon atom was developed by Wallace [281 ]. Slonczewski and Weiss [29]. 
A detailed single group analysis was given by Bassani and Parravicini [30]. Group theory 
has also been used [3l| to analyze the Raman scattering and electron-phonon interaction. 
Trigonal band structure was analyzed in terms of graphene double group [32.,] . Here, we 
ignore z plane reflection and consider the Cq v point single group, and the little group at K± 
as C% v single group. The symmetry group elements for C% v at Dirac points are E, C^,C 3 , 
<Jdi, where % = 1,2,3. The symmetry operation a^i leaves the Dirac points unchanged but 
interchanges A by B atoms in real space. a vi acts the opposite way. 

If we use to indicate the spinor reps for a rotation group, then for C% v , Ti ® D? = T 4 . 
We also have the relations r 5 <g> T 4 = r 3 , r 6 <g> T 4 = r 3 , and r 4 <g> T 4 = Ti + T 2 + r 3 . 
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TABLE I: Character table for double group C% v for K± point of graphene. T^s^ are extra repre- 
sentations (reps) and the spinor rep is T±. 



Without spin and when the Fermi level Ep crosses K± we have a T3 reps. Following Ref. 



29j we can show the Bloch sum for the p z orbitals on A and B sublattices are the basis of 
the representation, i.e. <p^(r) = -i= J2r x e ~* k '( Ri+TQ! )0"(r — Rj — r a ), where a = A or B, 
(f)f representing p z orbital at Rj + T a site. When spin is taken into account and spin orbit 
interaction (SOI) is ignored, the degeneracy at K± is 2 x 2 = 4. When SOI is considered, 
the reps of double group are split into r 3 <g> = T 4 + T 5 + r 6 . The reps T 4 , T 5 , and T 6 are 
degenerate in absence of SOI at K±. 

a) For ADA case the system has a C% v symmetry (different from C^ v , the character table is 
shown by table 1). The planes reflection symmetry in C^ v converts the point K + into K_ 
and vice versa. Thus, the symmetry C$ v mixes the states from the two Dirac points. In fact, 
this group is a subgroup of Cq v and the symmetry reflection planes do show up in Cq v and 
map Dirac points onto each others. Thus, only one flavor couples with the impurity. The 
Hamiltonian in this case shows that a basis belonging to Ti representation of C^ v can be 
constructed from the E rep of C$ v , which leads to an invariant hybridization Hamiltonian. 



If we consider explicitly a 3d magnetic impurity in this case, for example Mn (as realized 



in ADA case 



33|). the nonvanishing coupling is only present between the d z % orbital of the 



impurity and graphene since the d z 2 orbital transforms according to 1^ rep. 
b) For ADC case the symmetry group is Cq v where C 2 , and o vi show up in this group 

except the symmetry operations in C 3v . In this group, C 2 , Cq transform A sublattice into 
B sublattice, K + into K- and vice versa. a V i keep the sublattice unchanged but exchange 
K + and K_. Therefore, we have to construct the basis for Cq v by using the basis for the two 
Dirac points in C% v group. From Table 2, E rep in C% v corresponds to the E% and E2 reps in 
C Gv . Therefore we can use the Bloch sums (p^fka m the & re P * n ^° construct the reps 
in Cq v . Explicitly we consider a 3d magnetic element, and ignore the crystal electric field 
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TABLE II: Basis table of C§ v and C% v and their compatibility relations. R x = yp z — zp y is the 
angular momentum component around axis x, the others are obtained by cyclic permutation. 



splitting so that the energy levels for different reps are still degenerate or quasi-degenerate. 
The wave function should be constructed according to the symmetry group of Cq v . The 
reason is that the nonvanishing hybridization should be invariant under the operations in 
Cq v which is obtained as Tj £g> Tj by group theory. From Table 2, we infer d z % belongs to Ti, 
(d xz ,d yz ) to r 5 (_E' 1 ), and (d x 2_ y 2,d xy ) to T 6 (E 2 ). These 3d orbitals are well defined. Having 
specified the states of the impurity, we consider the basis of Cq v reps by combining the basis 



S 



of E in C 3v . We find 
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+ ^l_ka,0£_k«x 



+ ^0l +k J T , (6) 



where a; = e l2n ^ 3 . Therefore, the hybridization Hamiltonian for 3d magnetic impurity is 



#hyb = ( V ^ C lkafaa + h.C. 



(7) 



okcr 



where a = r^Ts,^ reps, and the latter two are 2D reps. The creation and annihilation 
operators are obtained by defining the field operators based on the respective basis of the 
corresponding rep in the standard way This situation resembles the Coqblin - Schrieffer 
model 34|, however the Hamiltonian now is written in terms of the irreducible reps of the 



systems. 

Kondo temperature under large U.- We calculated the Kondo temperature for the Anderson 
model in ADA case using the equation of motion method [351] that gives consistent results for 



the susceptibility with those derived by other methods at high temperature 



and delivers 



the correct Kondo temperature [37|. The Kondo temperature for large U is proportional 



to exp 



am ■ 



where A 



D 



t, D is the energy cutoff. If e < n and // — > 0, the Kondo 



temperature tends to zero. The impurity is then decoupled from graphene. If e < \x and 
with increasing the Kondo temperature increases exponentially. Thus, a small gate 



voltage may change the Kondo temperature dramatically 



9]. 



CONCLUDING REMARKS 



Summarizing, we investigated the Kondo model for a magnetic impurity in graphene in 
an adatom configuration. Starting from the tight-binding Anderson model we applied a 
partial wave method to expose the coupling flavors. The Kondo effect in graphene behaves 
as follows: i) For ADA case, only one flavor couples to the impurity, ii) For the ADC case, 
the number of channels is determined by the multiplets of the impurity as well. Multi-flavor 
character shows up in this case. Hi) Based on the irreducible basis of the reps of the 
system, A Coqblin-Schrieffer-like hybridization model is derived for a 3d magnetic atom. 
d z 2, (d xz ,d yz ), and (d x 2_ y 2,d xy ) couple with the r 1; T 5 (Ei), and Tq(E 2 ) reps respectively. 
The bases for these reps of graphene are derived explicitly, iiii) The degeneracy of the 
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two Dirac points only leads to a higher T K for ADA case. Applying a gate voltage the 
exchange coefficients are non-zero due to a finite DOS at E F and the Kondo effect should 
be observable. 

We thank C.-L. Jia, M.A.N. Araujo and A. Rosch for useful discussions. The work is 
supported by the excellence cluster " Nanostructured Materials" of the state Saxony- Anhalt. 
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